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ABSTRACT 


We find that a non-axisymmetric bow shock model, with an appropriate 
choice of parameters, could £t the line splitting velocity field of the cone-like Ha 
nebula in NGC 4945 better than a canonical cone model. Meanwhile, the bow 
shock model could also reproduce the morphology of the Ha nebula. The bow 
shock results from the interaction of the galactic superwinds with a giant HII 
region. It is implied that the starburst ring or disk around the galactic nucleus 
should be generating strong winds, and the bright Ha knot northwest of the 
nucleus be suffering an anisotropic mass loss process. 

Subject headings: galaxies: starburst — galaxies: individual: 

NGC 4945 — galaxies: kinematics and dynamics — shock waves 



1. Introduction 


It is now believed that the galactic superwinds may arise from the Seyfert, AGN, and 
starburst galaxies, and are able to create a giant bipolar wind cavity toward the galactic 
halo. The complicate phenomena in the cavity region are closely related to the nature of 
the superwinds and the galaxy itself. The edge-on spiral galaxy NGC 4945 has a wind 
cavity, and the observation of the split lines has been made (Heckman, Armus, & Miley 
1990, hereafter HAM). The P.A.= 45° line splitting region is ~ 500pc across, and the 
P.A.= 135° split lines are only found to begin from a distance ^70 pc to the galactic 
nucleus. Recently the conic wind cavity is displayed clearly in optical and near infrared 
wavelengths (Moorwood & Oliva 1994; Moorwood et al. 1996, hereafter MVKMO). The Ha 
image shows a conic nebula of an opening angle 78° inside an infrared fan-shaped structure. 
The nucleus is invisible in Ha, but is coincident with the position of the peaked emission of 
the K band and the Bry line. The spectrum of the head portion of the Ha nebula is typical 
of LINERS. 

It is suggested (MVKMO) that the cone-like Ha nebula can be explained with Suchkov 
et al.^s (1994, hereafter SBHL) theory on superwind bubble, in which the biconic cavity 
wall begins from a nozzle, and then stretches out like a cone. It could, however, be found 
difficult to fit the line splitting velocity field data with such a cone. For the position of the 
P.A.= 45° slit on the nebula, the cavity wall could be represented, as an example, by a 
cone described by HAM in concept, in which the materials flow out (at the speed Vq) along 
a biconic surface centered on the nucleus, with the cone’s symmetry axis along the minor 
axis of the galaxy. The origin of the coordinates is set at the central nucleus, and the axis 
is set along the minor axis, with the positive position toward the northwest. The aspect 
angle, namely the inclination of the line of sight to the galactic plane, is denoted as 0. If 
a distance d = 4.6 Mpc is adopted for NGG 4945 (HAM), and the ordinates of the cutting 
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point of the axis by the slit for the P.A.= 45° line splitting is written as (0, 0, Zq), then the 
projected distance of the cutting point to the nucleus is Zq coscj), measured 300 pc. With the 
parameters suggested by HAM, 0 = 33° and Vq = 525kms“^, the htting result is plotted 
as dotted lines in Fig.l. The curve for P.A.= 45° is of a feature something like a polished 
triangle, with a peaked top and a flattened bottom. This typical tendency will not be 
changed appreciably, if different values of parameters are adopted, and even if an expansion 
velocity component perpendicular to the cone surface is allowed. Though, the scenario that 
materials flow along a cone-like surface is constructive and indicative of another mechanism 
proposed as follows. 

On the other hand, as seen in the Ha image of the nuclear region of NGC 4945 
(Fig.Id & Fig.3a in MVKMO), there is not a visual “nozzle” as described by SBHL for 
this cone-like nebula, and the apex of this nebula is not located at the central nucleus, but 
rather nearly on the linking line between the nucleus and the bright Ha knot at the position 
“C” denoted there. The head portion of the nebula surrounds the bright Ha knot. This 
clear-cut cometary Ha nebula with a noticeable position of the Ha knot closely behind the 
nebular apex seems more likely to be a bow shock structure than a superwind cone. We 
hnd that not only the presence of a bow shock is possible because of the interaction of the 
galactic superwind with the Ha knot, but also the bow shock model could reproduce the 
morphology of the Ha nebula and £t the line splitting velocity held well. 


2. The bow shock model 

As we can see from the Ha image of the nuclear region of NGC 4945 (MVKMO), 
the bright Ha knot at “C” (denoted by MVKMO) has a nucleus-halo structure. And as 
HAM pointed out, the bright Ha knot at “C” is evidently a giant HH region, as well 
as other knots nearby. Based on the observation of the HH regions and the accordingly 
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established Champagne model for them (Yorke 1986; Osterbrock 1989), the HII region with 
a nucleus-halo structure should undergo a star formation activity {e.g., Ye 1992), and the 
stellar winds from the O stars in this type of HII region is undoubtedly important (see, e.g., 
Osterbrock 1989). We would assume that the kinetic energies of the winds of the individual 
O stars could be released out easily in spite of the collisions among themselves. Shock 
heating related to the ongoing starburst would be a dominant excitation mechanism for this 
region (see, e.g., Dopita & Sutherland 1995; Cozial 1996; Ho, Fillippenko, & Sargent 1993), 
and could be responsible for the LINER-like spectrum of region “C” obtained (MVKMO). 
According to Moorwood & Oliva (1994), the superwind of NGC 4945 is driven by starburst 
supernovae and/or AGN. The collision of the galactic superwind with the stellar winds 
from the HII region at “C” would inevitably produce a bow shock {e.g., Huang & Weigert 
1982; Bandiera 1993; etc.), and hence could account for the cone-like Ha nebula. Ha 
cometary nebulae have been very common in the stellar level and attributed to bow shocks, 
e.g. PSR 1957-1-20 (Kulkarni & Hester 1988), NGG 2359 (Schneps & Wright 1980; Ghen, 
Wang, & Qu 1995), etc., and the nebula in NGG 4945 seems very likely to be a bow shock 
phenomenon in the galactic level. In the following we would simulate the line splitting 
velocity held and the morphology of this nebula using a bow shock model. 

From MVKMO’s image, the projected distances of the apex of the Ha nebula to 
the knot at “G” and the nucleus measure = 53 pc and Za = 66 pc, respectively. If the 
superwinds arise purely from the nucleus, the ratio of these two distances ~ 0.8 would imply 
a much bigger opening angle (~ 160°) for the bow shock (Ganto, Raga, & Wilkin 1996). It 
should not be the real case, however. The Bry and L band images reveal a ~ 400 pc size 
starburst ring or disk on the galactic plane (Moorwood & Oliva 1994; MVKMO). Since this 
size is comparable with that of the line splitting region of the Ha nebula, the ring/disk can 
contribute a component of starburst driven winds, and then it could be expected to alter 
the whole effect of the superwinds. Such a wind complex could constrict the bow shock to a 
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smaller opening angle than the pure nuclear wind, though its momentum flux distribution 
is complicated. It is assumed for simplicity that the sum of the two components of winds 
can be approximated as a wind with an average velocity vi which acts as if it is generated 
from a convergence point Zi (< 0) far behind of the nucleus. 

As can be seen carefully in the Hq;+[NII] image, the head portion of the nebula is 
not axisymmetric with respect to the axis linking the nucleus and the knot at “C”: the 
nebular surface at the east is close to the knot and even a slight dent can be discerned, 
while the surface at the south is protrusive. This morphological dent could be due to the 
relatively low momentum flux of the giant HII region’s wind on its equatorial plane. When 
its equatorial plane is tilted, the asymmetry appears. Therefore, for the giant HII region, 
we would use a bipolar momentum flux distribution for the wind, which has been widely 
used in the stellar level. Furthermore, as the equatorial plane is tilted to the left, the left- 
and right-hand sides of the bow shock shell would swell up, and relatively the back- and 
front-sides would shrink inward (cf. the morphology simulation below). Consequently, the 
polished triangle configuration in the fitted velocity held for the axisymmetric cone would 
be changed, and especially the top part would be expected not to be peaked due to the 
flattened back side of the shell. An algorithm for calculating the asymmetric bow shock has 
been developed by Bandiera (1993), and would be adopted in our simulation as follows. 

The subscript 1 would refer to the quantities of the superwind complex, and 2 to the 
quantities of the wind from the Ha knot. The momentum loss per steradian from the Ha 
knot is presumably described by (Bandiera 1993) 

F2 = h;[l + Ay(l-3cos2x)/2]. 

Here is the average momentum loss per steradian. Ay (at a range —2 < Ay < 1) 
represents the amount of anisotropy in the wind momentum flux; when Ay < 0, 1^2 tends 
to be a bipolar distribution. And y is the angle between a given direction and the wind’s 
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axis of symmetry. The direction of the symmetry axis is represented by the azimuthal and 
polar angles (6^*,.^*) with reference to the axis. 

For the sake of making the freedom degree of the used parameters the least, the 
superwind complex is regarded as isotropic in the nebular region, and the velocity V 2 of 
the wind from the giant HII region as independent of directions hrst, which is presumably 
typihed by the FWHM ~ 600kms“^ of the broad feature in the spectrum of region “C” 
(MVKMO). Now the free quantities used to £t the data of the line splitting and the nebular 
morphology are Zi, Ui, 9 ^, Ay, and 0. Each of these quantities could sensitively affect 
some certain features in the outcome. Basically, zi affects the spatial range of the P.A.= 45° 
split lines; vi the blueshift and redshift values for the split lines; the projected position 
for the slight dented distortion in the nebular head portion; 6^* the projected position for 
the slight distortion in the head and the inclination of the closed curve for the line splitting 
velocity held. Ay the distortion degree of the head and the shape of the closed curve, and 
0 the height of the htting curves for the shifted velocities. Therefore there is actually very 
little freedom in the determination of the parameters, though the determination may not 
be unique. 

To calculate the blue- and redshift of the velocities, we project the streaming velocity 
along the bow shock surface to the line of sight. By repeated tests and comparisons, some 
parameters could be decided first: Zi — Za = —Ibrp, 6** = 180°, = 60°, and 0 = 14°, 

and then a solid and dashed lines are plotted in Fig.l for vi = 1400kms“^. Ay = —1.5 
and Vi = 1200 kms“^. Ay = —2, respectively. These two curves could basically trace all 
the points of data for the split lines, and are much better than the curve based on the 
cone model. A line splitting region ~ 500 pc at 2:0 cos 0 = 300 pc for the P.A.= 45° data 
is accordingly reproduced. The htting aspect angle 0 ~ 14° also seems reasonable for 
an almost edge-on galaxy (Moorwood & Oliva 1994). And it is in better agreement with 



that observed for the galaxy (12°) by Nakai (1989), than that estimated by HAM (33°). 
With the ham’s estimation of the mass loss rate of the nuclear region of NGC 4945, 
1.6M0yr“^, the htting velocity vi ~ 1400—1200kms“^ would lead to a momentum flux 
of ~ 1.4—1.2 X 10^^ dynes. This is in a very good agreement with HAM’s estimation of 
the momentum flux ~1.3, 2, or 5x10^^ dyne s“^ in different ways. Major shortcomings of 
the above simulation may be that 1) there is a small bump in both the htting curves for 
P.A.= 135°, and 2) the two curves for P.A.= 45° could not reproduce the trend of a little bit 
larger splitting of the southwestern half. They could be understood if we consider that the 
real environment of the winds are more sophisticated than described in our relatively simple 
model. For example, the small bumps in the htting curves for P.A.= 135° could possibly 
be accounted for with the ignorance of the additional squeezing ehect of the starburst 
ring’s wind to the head portion of the bow shock in the calculation, because this additional 
ehect could make the head portion more acute than the present simulation. The L(3.5 /rm) 
continuum image (MVKMO) and the bumps in the H 2 rotation curves seem to be more 
suggestive of a ring than a disk (Moorwood & Oliva 1994). And the actual bow shock tip 
might be closer to the nucleus than the observed apex of the Ha emission, and be heavily 
extincted; thus the bump could be expected to be alleviated. The latter situation would 
suggest that the wind from the nucleus itself is not as strong as the wind from the extended 
part (he., the ring/disk). Considering the asymmetry in the galactic wind complex could 
possibly remove the bump and £t the data better, too. As for the second point, it could 
be improved by introducing an anisotropic velocity distribution for the wind from the HIT 
region: 

V 2 = V 2 a{^ - cosy), 

where V2 is assumed to have the same symmetric axis as Ay, V2a is the average velocity of 
the wind (~ 600kms“^), and represents the amount of anisotropy in the wind velocity. 
While maintaining the parameters adopted above, we take A^ = 0.5. The resulting split 
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velocities are plotted in Fig.2. The curves for P.A.= 45° fit the data points and the 
mentioned trend pretty well. 

The morphology of the bow shock is decided by the momentum flux, and is independent 
of the wind velocities. Taking a set of parameters used above Zi — Za = —16ro, 6^ = 180°, 
= 60°, Ay = —1.5, and (j) = 14°, for instance, the morphology of the bow shock is 
demonstrated in Fig.3, which reproduces the shape of the Ha nebula very well. 

In view of the above simulations, an asymmetric bow shock model could indeed £t 
both the line splitting data and the nebular morphology. In the simulation, the galactic 
gravitation has been ignored. In fact, one could hnd that the gravitation effect is not 
important in the head portion compared with the wind’s ram presure, and that the flow 
velocities (^500kms“^) in the trailing portion is higher than the escape velocity. The 
reasons why the stellar winds from the giant HII region could have anisotropic momentum 
flux and velocity distributions remain to be further studied. Here a preliminary suggestion 
may be that the HII region is flattened, with denser materials on its equatorial plane than 
at the two poles. 


3. Discussion 

The fan-shaped appearance of the hlaments emerging from the nucleus of NGC 4945 
discovered by Nakai (1989) is an evidence for the existence of a superwind activity. The 
opening angle of this fan structure (~ 150°) is much larger than that of the cone-like Ha 
nebula (~ 78°) only within which the split lines were detected (HAM). This fact implicates 
that the filament fan may have a common character with the infrared fan-shaped structures 
in Z, J bands, in which the extinction is lower (MVKMO). So the bow shock Ha nebula 
is formed within the superwind cavity and is different in essence from the fan-shaped 
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structures. The interior of the cone-like Ha nebula should be occupied by the wind from 
the HII region, rather than the galactic superwinds. 

It is intriguing that, in the Ha image of NGC 4945, there is another bright Ha knot 
near the denoted position “B” on the opposite side of the “C” knot, aligned symmetrically 
with respect to the nucleus. It is also worth noting that the extension of both [Fell] and 
H 2 ( 1 - 0 )S'( 1 ) line emission are roughly symmetric on both sides of the galactic plane (as 
pointed out by Moorwood & Oliva 1994). These phenomena of symmetry might have 
something to do with the symmetric activity from the nucleus, most likely a bipolar galactic 
superwind. MVKMO suggested that the infrared fan structure is consistent with the 
superwind bubble model developed by SBHL. Indeed, the two bright HII regions discussed 
above could naturally be born as a result of the bipolar superwind. According to the SBHL 
model, especially to the B2 model therein, the superwind entrains the disk matter to the 
halo. At very early stage, two dense gaseous blobs, or two clouds, are formed symmetrically 
about the galactic nucleus. Thus two giant HII regions, as we see in the nuclear region of 
NGC 4945, might thereafter be hatched in these two dense clouds due to the star formation 
activities, giving rise to the stellar winds. 


4. Conclusion 

We have found it difficult to £t the line splitting velocity held of the conic Ha nebula 
using a canonical cone; with an appropriate choice of parameters, however, it could 
essentially be reproduced by a non-axisymmetric bow shock. The morphology of the bow 
shock simulated is also very similar to the Ha nebula. It is implied that the the starburst 
ring or disk around the galactic nucleus should be blowing strong winds. The bow shock 
originates from the collision of the superwind complex from the nucleus and the starburst 
ring/disk with the northwestern bright Ha knot, namely with the wind from the giant 
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HII region at “C”. The superwind complex could be effectively approximated as a wind at 
an average velocity ~ 1400 km s“^ from a convergence point far behind the nucleus. And 
the wind from the giant HII region might have anisotropic momentum flux and velocity 
distributions. The best-fit average velocity of the superwind and the aspect angle are in 
good agreements with the observational data. 

We wish to thank Rino Bandiera and Tim Heckman for careful reading and critical 
comments for the manuscripts which help to improve our analyses. We are indebted to R. B. 
also for providing his elaborate codes for computing the non-axisymmetric bow shock. This 
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Analysis of the Department of Astronomy, Nanjing University. This work is supported by a 
grant from the NSF of China, a grant from the State Education Commission of China for 
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Fig. 1. — The simulation of the line splitting velocity field, the top hgure for P.A.= 45° and 
the bottom hgure for P.A.= 135°. The crossed points are adapted from RAM’s observational 
data. The dotted line is plotted according to the cone model with RAM’s parameters except 
that the apparent opening angle 9app = 78° has been translated to the real angle 9ri by 
projectional effect: 9ri = 2 arctan[tan(6'app/2) cos0]. The solid and dashed lines are obtained 
from the bow shock model for Zi = Za — IGr^, Vi = 1400 kms“^, 9^ = 180°, = 60°, 

Ay = —1.5, 0 = 14°, and zi = Za — 16ro, vi = 1200 kms“^, 9^ = 180°, = 60°, Ay = —2, 

0 = 14°, respectively. (See the text for the denotation of these quantities.) 

Fig. 2.— The simulation of the line splitting velocity held using the bow shock model with 
anisotropic distribution for the RR region’s wind velocity. The top hgure is for P.A.= 45° 
and the bottom hgure for P.A.= 135°. The parameters for the solid and dashed lines are the 
same as in Fig.l, except that A„ = 0.5 in addition for the both cases. (See the text for the 
denotation of the quantities.) 

Fig. 3.— The projected morphology of the bow shock for zi = Za — 16ro, 9^ = 180°, = 60°, 

Ay = —1.5, and 0 = 14°. (See the text for the denotation of these quantities.) 



